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Abstract 



m 

We study the correlation of the velocity dispersion of the coexisting molecules H 13 CN 
and H 13 CO + and the turbulent energy dissipation scale in the DR21(OH) star-forming 
region. The down-shift of the H 13 CO + spectrum relative to H 13 CN is consistent with 
the presence of ambipolar diffusion at dissipation length scales that helps the process of 



i-G turbulent energy dissipation, but at a different cut-off for ions compared to the neutrals. 

We use our observational data to calculate a turbulent ambipolar diffusion length scale 
L ~ 17 mpc and a strength of B pos ~ 1.7 mG for the plane of the sky component of 

r/3 , the magnetic field in DR21(OH). 
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" . ■ 1. Introduction 
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Interstellar magnetic fields and turbulence have been the most debated mechanisms of support 
against gravitational collapse in molecular clouds. Recent studies have evol ved towards acknowl- 



edging both mechanisms by studying turbulence in a magnetized gas (e.g., IShebali n et al. 



Goldreich fc Sridhar|[l995l ; bstriker et aDboOll : iBiskamplhocd : IfJasu et all2009l : lFalceta-Goncalves et al 



1983; 



2010J). It is, however, a challenge to combine the existing models for turbulence with observational 



data for the measurement of the interstellar magnetic fields to fully understand the dynamics in 
molecular clouds. There is yet no closed theory for magnetized turbulence due to its complicated and 
intermittent nature. On the other hand, interstellar magnetic fields are weak (~ few hundred /xG 
to a few mG) and difficult to measure with the very few observational techniques currently at hand 
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(e.g., ICrutcher et al.lll999l : iHoudel 12004 : iFakarone et all 12008 : iHildebrand et al.1 12009 : iHoude et al 
2009). It is therefore essential to innovate with more observing techniques to provide reliable data 
and test existing models. 

Turbulence is basically a dissipative process. In an incompressible turbulent fluid, energy is 
transferred from la rge eddies to smaller ones, until it is dissipated by viscosity at small enough scales 
(|Kolmogorovl Il94ll ). The dissipation of turb ulent energy has important effects on the ISM, such 
as initiating collapse in supercritical cores (jBasu et al.l [2009j) , modifying the chemical evolution 
of the gas molecules (IXie et al.l 1 19951 ) , and causing the decoupling of gas and grains in a cloud 
(JFalgarone fc Pugetlll995J ). 



Li fc Houde (120081) estab lished a new observational technique, expanding on the original work 



of Houde et al. (J2000al : l2000bl ) on the effect of magnetic fields on the spectral line widths of ion and 
neutral molecular species, to calculate the ambipolar diffusion length scale and the strength of the 
plane-of-the-sky component of the ambient magnetic field (B pos ) in molecular clouds. They studied 
the turbulent velocity dispersion spectra of a pair of coexistent ion and neutral molecular species 
(HCN and HCO + ) in M17, with the assumption that the turbulent energy dissipation process is 
associated with ambipolar diffusion. 



In this paper, we explain the method of iLi fe Houdd (120081 ) and our purpose to test their 
turbulent energy dissipation model using optically thin coexistent ion and neutral molecular species 
in §2, and present our observations in §3. The data analysis is explained in §4 and we end with a 
summary in §5. 



2. The Method of ILi & Houdd ([2008) 



Ostriker et al.l (|200ll ) introduced a numerical method to study the kinetic structure in a single- 



fluid turbulent cloud by measuring the distribution of the velocity dispersion against the size and 
mass of the simulated regions over which it was averaged. They showed how the lower envelope 
of the calculated velocity disp e rsions depicts the actual turbulent velocity spectrum of the gas in 
their simulations. ILi &: Houdd (120081 ) used this technique to study the velocity dispersions of HCN 
(J = 4 — >• 3) and HCO + ( J = 4 — > 3) spectral line profiles in M17 as a function of length scale. 
They noticed that the lower envelope of the HCO + data was consistently down-shifted from that 
of HCN by a constant amount, independent of scale, and interpreted this result as the signature 
of turbulent ambipola r diffusion ac t ing at smaller spatial scales than probed by their observations. 
Based on these results ILi fc Houdd (|2008l ) suggested a scenario where ambipolar diffusion sets in at 
small enough length scales where neutrals decouple from the flux-frozen ions and the magnetic field, 
and the friction betwee n the ions and the drifting neutrals causes the dissipation of turbulent energy. 
Tilley fc Balsaral (12010 ) tested this model by simulating a two-fluid (ion and neutral) turbulence with 
a range of different ionization fractions. Their results showed that at length scales comparable to the 
ambipolar diffusion scale certain MHD waves, including Alfven waves, are strongly damped, and this 



damping causes the velocity spectrum of ions to drop below that of the neutrals. They simulated ion 
and neutral spectral line profiles for ionization fractions ranging from 10~ 2 to 10~ 6 and plotted their 
corresponding velo city dispersions against different beam sizes, producing a result very similar to 



Li &: Houdd (120081). Thus they were successful in numerically duplicating the observational results 



of lLi fe Houdd (J2008J ) . including their determination of the turbulent ambipolar diffusion scale. 



Assuming the gas to be well coupled to the field lin e s at l arge scales and that ion-neutral 
decoupling occurs at dissipation length scales, iLi &; Houdd (|2008l ) obtained an expression for B pos 



by setting the effective magnetic Reynolds number R m ~ V n L//3 to 1 at the decoupling length scale. 
Here (3 = B 2 / 'Aimi/Ai/i is the effective magnetic diffusivity with n, the ion density, Vi the collision 
rate of an ion with the neutrals and \i the mean reduced mass characterizing such collisions. The 
aforementioned expression for B pos is 

/ jl \ 1/2 / y , \ 1/2 n l/2 

Bpos = [oJ^Tc ) ^lkms-ij (l06cm-3j (l0=0 mG ' (1) 

where L is the ion-neutral ambipolar diffusion decoupling length scale, V n the neutral velocity 
dispersion at L , n n the neutral volume density, and Xe is the fractional ionization. It should be 
noted that Equation (p} provides only an approximate value for B pos as it assumes that ambipolar 
diffusion sets in exactly when R m = 1, which is most likely too simple an approximation. The 
decoupling scale and the velocity dispersion V n are determined by fitting the lower envelopes of the 
square of the ion and neutral velocity dispersion data (a 2 and a 2 , respectively) to Kolmogorov-type 
power laws. If the coexisting ions and neutral species are well coupled at large scales, their turbulent 
power spectra will have the same spectral index over the inertial range (i.e., at scales larger than 
the dissipation scale and smaller than the injection scale). Moreover, the cut-off range for the ion 
and neutral energy dissipation do not coincide, and assuming the ion energy spectrum to drop 
significantly at the ambipolar diffusion scale, the difference between a 2 and a 2 is the difference of 
the integral of the ion and neutral energy spectra between the ambipolar diffusion scale and neutral 
vi scosity dissipation scale (this difference, a (see Equation [2] below), is the shaded area in Figure 3 



of ILi &: Houdd 120081 ) . With all the above assumptions, the a 2 and <r 2 data are thus fitted to the 



following expressions 



a 2 (L) = bL n + a (2) 

a 2 n (L) = bL n . 

The observational determination of velocity dispersions from spectral line profiles is a measure- 
ment of a velocity field that is integrated along the line of sight within a column of gas of volume 
I 2 x L, where I is the beam size and L is the depth through the cloud. However, a given velocity field 
is associated with turbulent substructures that span volumes different than I 2 x L. It should there- 
fore be examined whether the observational value of a 2 at a given beam size / truly represents the 



actual velocity dispersion at that length scale. iFalceta-Goncalves et al.l (J2010I ) investigated this issue 
with spectral line profile synthesis and three-dimensional MHD numerical simulations. They showed 
that although the minimum observed velocity dispersions scale as the actual a 2 , the coincidence 
of individual values at a given scale depends on I and the sonic Mach number of the system. The 
observed dispersion is generally a good approximation of the actual dispersion, although the latter 
is slightly underestimated for supersonic turbulence and overestimated for a subsonic regime. They 
also showed that the power law functions fitted to the lower envelopes and actual three-dimensional 
velocity dispersions from the simulations correspond well. The dispersions obtained using the lower 
envelope therefore appear to be reliable estimates for the actual velocity dispersions, as long as a 
large enough number of lines of sight are considered. 



Li et al.l ( 20101 ) studied the differences in the velocity dispersion for th e spectral l ines of H 13 CN 
(J = 1 -». 0) and H 13 CO+ (J = 1 -)• 0) in DR21(OH) using the data of lLai et all J2003aJ ). They 
performed a similar analysis for HCN (J = 3 — > 2) and HCO + (J = 3 — ■> 2) data they obtained 
with the Submillimeter Arra y in NGC 20 2 4, an d also HCN ( J = 4 - » 3) and HCO+ (J = 4^3) 
lines previously observed by iHoude et al.l (J2002J) in M17 and used in iLi fe Houdd (120081 ) . In each 
source, the velocity dispersion comparisons were performed for two locations with measured B pos 
line densities. They found that the ion-neutral velocity difference is almost constant in M17, which 
has a uniform density of magnetic field lines, while in NGC 2024 and DR21(OH) the ion- neutral 
velocity dispersion difference correlate wel l with the magnet ic field line density. This behavior is 
predicted and consistent with the model of lLi &; Houdd (J2008I ). They also considered other possible 
parameters that could cause differences between the velocity dispersion of the two tracers, including 
optical depth, outflow enhancement, and presence of hyperfine structures in the spectral lines, but 
only turbulent ambipolar diffusion could consistently explain their observations. 

The spectra of an HCN and HCO + pair, however, usually appear to be optically thick and 
ar e subject to s aturat ion. The aim of this work is to test the turbulent energy dissipation model 



of ILi fc Houdd ( 20081 ) with optically thin isotopologues i.e., H 13 CN and H 13 CO + , to investigate 



whether the constant down shift of the velocity dispersion spectrum of the ion from that of the neu- 
tral is consistently reproducible with different species. Accordingly, the lower optical depths of the 
aforementioned molecules will allow us to pote ntially probe magnetic fields at greater depths in star- 
forming regions, where they may be stronger (ILi fc Houddl2008l : lBasull2000l : ICrutcherlll999l ). Also, 
just like their main isotopologues they have very similar molecular masses and the critical densities 
associated with their J = 4 — » 3 transition are approximately 10 cm -3 and 10 6 cm -3 , respectively. 
In particular, their line profiles will not be saturated or self absorbed, which is likely to be the case 
with HCN and HCO + , and therefore the calculated line widths will be better representatives of the 
true velocity dispersions characterizing the medium under study. 



3. Observations 

We obtained position-switch grid maps of H 13 CN(J = 4 -> 3) at 345.336 GHz and H 13 CO+( J = 
4 -> 3) at 346.998 GHz in weather grade 3 (0.08 < r 22 5 < 0.12) in DR21(OH) during the months 
of August to November 2008 at the James Clerk Maxwell Telescope (JCMT ) on the summit o f 



Mauna Kea, Hawaii. The observations were performed using the HARP-B (jBuckle et al.l 120071 ) 
single sideband receiver and the ACSIS correlator configured for a bandwidth of 500 MHz and 
channel resolution of 61 kHz. HARP-B is a 4 x 4 element heterodyne focal plane array that uses 
SIS detectors. The 16 detectors have receiver temperatures of 94 — 165 K and are separated by 30" 
with a foot print of 2'. To fully sample the 15" JCMT beam, we mapped with sample offsets of 15" 
around the reference position (a = 20 39 m 01 s and 5 = 42°22 37.7 , J2000) of the source. Since the 
difference between the frequencies of H 13 CN(J = 4 -> 3) and H 13 CO+(J = 4 -)• 3) is 1.66 GHz, we 
made use of the special configuration capability of ACSIS to observe the two spectra simultaneously 
within the 1.8 GHz passband of the receiver. We used the previousl y published main be am efficiency 



of ~ 62 % determined with observations of Saturn at 345 GHz by iBuckle et al.l (J2007I ). 



4. Results 

Figured] shows the intensity maps of H 13 CN ( J = 4 — > 3) and H 13 CO + ( J = 4 — > 3), integrated 
over a velocity range of —20 kms -1 to 20 kms -1 in DR21(OH). The high critical density of the 
J = 4 — > 3 transition results in detection of radiation from high density regions only, as can be seen 
from the intensity contours congregating in an area of ~ 40" x 40" around the reference position 
of the source. We performed the spectral line analysis with the Starling and Gildaaj software 
packages. The spectra of neighboring regions were averaged across the maps to simulate beam 
sizes 15", 30", 45" and 60" at different lines of sight in the observed region. We then fitted multi- 
Gaussian profiles to each spectrum and the corresponding velocity dispersion a was then calculated 
using these fits. The plot of the square of the velocity dispersions for which the dispersion is at 
least three times larger than its corresponding uncertainty is shown in the top panel of Figure [2] 
as a function of length scale, with the H 13 CN data plotted in black and H 13 CO + in red. We did 
not include the error bars in this plot for a clearer representation of the data points. The H 13 CN 
data corresponding to the largest beam size are not shown, as they do not fit in the limits of the 
plot; we had low statistics at the 60" scale and therefore were not able to produce quantitative 
data points for either of the species at this scale. The presence of data points at each length scale 
above the lower envelopes can be associated with the interception of the corresponding line of sight 
with several turbulent substructures of different velocities that brings an increase in the observed 
velocity dispersions. Therefore, the minimum values for a 2 may be attributed to the line of sight that 
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intercepts the fewest number of such turbulent sub-structures (jFalceta-Goncalves et al.ll2010l ). All 
the minimum a 2 values we obtained for both species arise from one location in DR21(0H), circled 
in red in Figure [H thus further confirming the coexistence of H 13 CN and H 13 CO + and showing 
the specific region along the line of sight that contains the least number of turbulent substructures. 
The spectra for the different beam sizes at this location are displayed in Figure [3l with the H 13 CO + 
profiles plotted in red and H 13 CN in black and scaled to the temperature of H 13 CO + for a clearer 
comparison of their line widths. 

The lower envelopes of a 2 for H 13 CN and H 13 CO + are plotted with their corresponding uncer- 
tainties in the bottom panel of Figure [21 with Kolmogorov-type power law fits given by Equation 
([2|). As mentioned above, we did not have sufficient statistics to find the true lower envelope at 
the largest scale (i.e., 60") for both spectra and thus used the a 2 data at the three smallest beam 
sizes for the power law fits. To enforce the same power law fit for both species, we first fitted the 
difference between the square of the velocity dispersion data of H 13 CN and H 13 CO + to a constant 
function to obtain a, and then fitted the sum of the two data sets to a power law of the form 
2bL n + a to obtain b and n. The fit results provide a = -0.45 ± 0.01 km 2 s -2 , b = 0.49 ± 0.03 
km 2 s -2 arcsecond - ™ and n = 0.36 ± 0.02. Altho u gh ou r fit is not as equally good for the two 



species at all length scales as it was for iLi fe Houdd (|2008l ) , it is still very satisfactory and perfectly 
adequate for our study. Assuming a significant drop in the ion velocity spectrum at the dissipation 
scale, the velocity dispersion of ions w ill be approximately equal to V n at L i.e., V n (L ) ~ <r 2 (L ). 
We can then write (JLi fc Houdd bood ) 



V^{L') ~ 0.37bnL' n 
~ a + bL 



(3) 



and therefore 



L 



-a/[b(l - 0.37n)]. 



(4) 



or 



Taking the distance from DR21(OH) to be 3 kpc (JGenzel fc Downeslll977l ). we obtain V 
i 17 mpc and V n ~ 0.26 km s — 1 . 



1.2" 



We used the RADEX software package (IVan der Tak et al.l 120071 ) to calculate the gas volume 
density and optical depths associated with the J = 4 — > 3 transition of H 13 CN and H 13 CO + . 
Taking the values for the line widths and line temperatures of the spectral lines observed at the 
reference position of the source, and the column den sities and excitation temperatures of the two 



species (~ 10 13 cm 2 and ~ 13 K, respectively) from iHezareh et al.l ( 20081 ) . we get n n ~ 10 6 cm 3 
and assuming the observed core to be as deep as it is wide, i.e., 40" , we obtain A(H2) ~ 2 x 10 24 
cm" 3 . Additionally, RADEX computes optical depths of 0.17 and 0.21 for H 13 CN and H^ 3 CO+ lines, 
respectively. This optical depth comparison agrees with the discussion by lLi et al.l ()2010l ) mentioned 



earlier. We also use Xe = 3.2 x 10 8 for the ionization fraction as obtained bv lHezareh et al.l ( 20081 ) 



for DR21(OH). Applying the parameters above in Equation ([I]), we determine a value of B pos ~ 1.7 
mG. This value can be inaccurate by a factor of a few due to the uncertainties in n n and Xe, and the 
fact that Equation (JTJ) is derived with the aforementioned assumption of determining the turbulent 
ambipolar diffusion length scale by setting R m = 1. 



Our results are in general consistent with the calculations of iLi fc Houdd (|2008l ) although 
the measurements were done on different sources and with different pairs of molecular species. 
An exception is pe rhaps the measured turbulent ambipolar dissipation scale for Ml 7 reported by 
Li &: Houdd (|2008l ) was 1.8 mpc, almost an order of magnitude smaller than what we measure for 
DR21(OH). This, however, may simply be the result of a stronger magnetic field for our source. 
Indeed, using the same fractional ionization of 3.2 x 10 -8 for M17 yields B pos ~ 0.57 mG for that 
source, approximately three times weaker than for DR21(OH). A relatively weaker field strength for 
Ml 7 than for DR21(OH) is also corroborated with existing CN (N = 1 — >■ 0) Zeeman measurements 



(Crutcher et al. 



that of 



Lai et al. 



1999 



Falgarone et al.l 12008). Furthermore, our result for B pos is consistent with 



(2003b), wh o obta ined a polarization map of dust emission in DR21(OH) and used 



the IChandrasekhar fc Fermil (J1953I ) method to calculate B pos — 1 mG. 



We can now c ombin e our value for -Bp OS with the Zeeman measurement of B\ os ~ 0.4 mG 



by ICrutcher et al.1 (|l999) in DR21(OH) and obtain a total strength of B ~ 1.8 mG at a volume 
density n(R2) — 10 6 cm" 3 . We can also use this result to calculat e a mass to flux ratio for the 
observed source. Taking M/$ B = 1.0 x 10~ 20 -/V(H 2 )/|.B| cm 2 /xG (JMouschovias k. Spitzerlll976l ) 
with N(R 2 ) ~ 2 x 10 24 cm" 2 , we obtain M/<S> B ^ 10, which implies that DR21(OH) exhibits 
a highly supercritical core. Although this number could be overestimated by the column density 
considering our assumption for the core geometry, it is consistent with the f act that DR21(OH ) 
has indeed re yealed signs of high-mass star formation. Detection of OH masers (JNorris et al.lll982l ). 
H2O mase rs (IGenzel &: Downed 119771 ). and high-velocity outflows from its two dense cores, MM1 
and MM2 (jLai et al.ll2003bl ). are indications of ongoing star formation in this source. 



5. Summary & Discussion 

We simultaneously observed and mapped H 13 CN (J = 4 — » 3) and H 13 C O + (J = 4 — > 3) in 
DR21(OH) to test the turbulent energy dissipation model of lLi fc Houdd (120081 ) . The down shift of 
the ion velocity dispersion spectrum from that of the neutral is readily explained by the presence of 
ambipolar diffusion that causes a steeper t urbulent power spe ctrum for the ions compared to that 
of the neutrals at dissipation length scales (jLi fc Houdd 12008). We used our observational data to 
calculate a turbulent ambipolar diffusion scale of L ~ 17 mpc and a plane of the sky magnetic field 
strength of B 



pos 



1.7 mG for this source. 



Our analysis shows that the difference between the velocity dispersion of the optically thin 
H 13 CN and H 13 CO + spectral lines at different length scales exhibits a trend similar to the results 



originally obtained bv lLi fc Houdd (J2008I ) with the optically thick main isotopolo gues. It therefore 
follows that t he ion line n a rrowi ng effect first discussed by iHoude et al.l (J2000al ) and the ensuing 
conclusions of iLi &: Houdd (|2008l ) cannot be due to a relative op tical depth ef f ect bet ween the ion 
and neutral species. This further corroborates prior re sults from IHoude et al.l (|2000bl ) . which first 
investigated this possibility, and more recently those of ILi et al.l (|2010l ). 



Although, as previously mentioned that optically thin spectral lines are preferable for this 
analysis, it should be noted that H 13 CN and H 13 CO + are less abundant than the main species and 
are more difficult to map. Indeed, we observed the aforementioned pair in three additional sources 
DR21(Main), NGC2264 and NGC2068 but DR21(OH) was the only source in which we obtained 
enough detections away from the core to perform our analysis. Moreover, observing HCN and 
HCO + in regions where they are optically thick in the line core should not be too detrimental to 
the line width calculations, as the wings of such spectral lines are still optically thin and suitable for 
this study. More precisely, since the wings are more heavily weighted in line width calculations, the 
(multi-) Gaussian profiles fitted to such sections, away from the line peaks can be used for obtaining 
the velocity dispersions in both species (|Li et al.ll2010l ). 



Finally, this technique has the advantage that it is not restricted to nearby sources, since the 
data analysis mainly involves line width calculations and does not require a spatial resolution high 
enough to resolve the ambipolar diffusion scale. 

The authors thank the referee for providing helpful comments, and also S. Basu and W.Dapp 
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Research Chair, Canada Foundation for Innovation, Ontario Innovation Trust, and Western's Aca- 
demic Development Fund programs. 
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Fig. 1. — Comparison of the intensity maps of H 13 CN and H 13 CO + ( J = 4 — >■ 3) integrated over a 
velocity range of -20 kms -1 to 20 kms -1 in DR21(OH). The contours span a range of 10% to 90% 
of the peak (9.03 Kkms" 1 for H 13 CN and 8.23 Kkms" 1 for H 13 CO+) by increments of 10%. The 
red circles display the location along the line of sight where the minimum values of a 2 are obtained 
observationally. 
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Fig. 2. — (Top) Plot of a 2 as a function of length scale in DR21(OH). The error bars for the 
data points are not shown to avoid clutter in the graph. The values of a 2 for H 13 CN at 60" fall 
above the plot margin. (Bottom) The lower envelopes of the H 13 CN and H 13 CO + data fitted for 
Kolmogorov-type power laws. 
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Fig. 3. — Comparison of the line widths of the J = 4 — > 3 transition of H 13 CO + (red) and 
H 13 CN (black) spectral lines corresponding to the minimum a 2 values. Starting from top left, going 
clockwise are the spectra at beam sizes 15", 30", 45" and 60". In all the figures, the H 13 CN lines 
are scaled to the H 13 CO + line temperatures to provide a clearer comparison of the line widths. 



